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Abstract 



We present the first energy and angle resolved measurements of e'^'e^-pairs emitted 
from heavy nuclei {Z > 40) at rest by internal pair conversion (IPC) of transitions with 
energies of less than 2 MeV as well as recent theoretical results using the DWBA method, 
which takes full account of relativistic effects, magnetic substates and finite size of the 
nucleus. The 1.76 MeV EO transition in ^^Zr {^^Sr source) and the 1.77 MeV Ml transi- 
tion in ^°^Pb (^°^Bi source) have been investigated experimentally using the essentially 
improved setup at the double-ORANGE /3-spectrometer of GSI. The measurements prove 
the capability of the setup to cleanly identify the IPC pairs in the presence of five orders 
of magnitude higher (3~ and 7 background from the same source and to yield essen- 
tially background-free sum spectra despite the large background. Using the ability of the 
ORANGE setup to directly determine the opening angle of the e+e^-pairs (Og+e-), the 
angular correlation of the emitted pairs was measured within the range covered experi- 
mentally (40° < ©e+e- < 180°). In the ^°Zr case the correlation could be deduced for a 
wide range of energy differences of the pairs (—530 keV < Ej\ < 530 keV). The ^°Zr 
results are in good agreement with recent theory. The angular correlation deduced for 
the Ml transition in ^°^Pb is in strong disagreement with theoretical predictions derived 
within the Born approximation and shows almost isotropic character. This is again in 
agreement with the new theoretical results. 

PACS numbers: 14.60.Cd, 23.20.En, 23.20.Ra, 29.25.Rm, 29.30.Aj 
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1 Introduction 



Although known since the 1930's as one of the elementary decay channels of exited nu- 
clear states with transition energies of more than 2meC^ = 1.022 Me V JI], 0], only limited 
theoretical and very limited experimental investigations of this process have been car- 
ried out in the following years. Theoretical results concerning the angular correlation of 
the e"''e~- pairs were mainly based on Born approximation and thus valid only for point- 
like \ow-Z nuclei and/or transition energies much higher than the threshold of 1.022 MeV 
0, ^, 1^, 1^, 1^, 1^. Only in the early 1990's more theoretical effort yielded results valid also 
for finite size nuclei and low transition energies (below twice the threshold) at least for 
electrical multipole transitions 0. Correct treatment of magnetic transitions was achieved 
only recently |T0[, and the results for angular correlations are presented here. 

Very few experimental results on IPC are published, and most of them on total IPC 
coefficients only, derived from the integral e"*" yield [|ll|. However, one very early experi- 
ment ]T^ , [T^ already measured the angular correlation for the 6 MeV El transition in ^^O, 
but without any energy measurement. This and other transitions of even higher energies 
in light nuclei have been measured recently fl^ both energy and angle resolved. The 
results were in reasonable agreement with Born approximation, as expected, since this 
approximation should be valid for these cases. However, e^e~-pairs emitted via IPC from 
heavy nuclei are of considerable interest, since IPC is the only known source of e^e~-pairs 
with a sharp sum energy. Narrow lines with sum energies {Es = E(.+ + E^.- ) in the range 
500 keV < Ey, < 800 keV have been found in heavy-ion experiments around the Coulomb 
barrier at GSI |T5|, |T6|, |1^, |18|, |T9|, The angular correlation of the emitted leptons 

is of particular interest, since the response of the various setups depends strongly on it. 
Prior to this work, no data have been available. Additionally, the most recent theoretical 
results shown below should be tested experimentally. 



2 Theory of internal pair conversion 

The angular correlation of electron-positron pairs emitted by internal pair conversion (mul- 
tipolarity L > 0) can be expressed in terms of Legendre polynomials 
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, ^ , r = - -TT^ 1 + > aiPAcosv)]. 1 

dEdcosT? 2 dE V ^ * 7 ^ ^ 

The differential pair conversion coefficient ^ and the anisotropy coefficients are 

calculated numerically. The calculation is performed using DWBA, i.e., the electron and 

positron wave functions are taken to be the exact scattering solutions of the Dirac equation 

for the Coulomb potential of an extended nucleus The differential pair conversion 

coefficient and the anisotropy coefficients depend therefore on the nuclear charge number, 

as well as on the nuclear transition energy and the multipolarity and parity of the nuclear 

transition. Thus, the calculation goes beyond the Born approximation, where electron and 



positron wave functions are considered as plane waves B. 41; is also tabulated in |T0 



Fig. 1 shows the angular correlation for Ml pair conversion of an lead nucleus {Z = 82) 

and a transition energy of 1.77 MeV. The energy difference Ea = Ef.+ — E^.- amounts 

to 150 keV, which corresponds to the experimental data presented in this work (see Fig. 

10 below). The calculation was performed assuming a point-like nucleus (dashed curve) 

as well as an extended nucleus (full curve). The point nucleus approximation is not very 

realistic for magnetic pair conversion in highly-charged nuclei since it overestimates the 

pair emission rate. Finally the dotted curve corresponds to the result obtained within the 

plane wave Born approximation [^], which displays the typical forward-peaked behaviour 

previously expected for the pair emission. For high nuclear charge numbers, the angular 

correlation deviates strongly from the Born approximation result. The pair emission occurs 

then nearly isotropic. This behaviour is elucidated by Fig. 2, which depicts the angular 

correlation, assuming a 1.77 MeV Ml transition with symmetric splitting of the transition 

energy on the electron and the positron and various nuclear charge numbers. 

Also in the case of electric monopole pair conversion (denoted as EO pair conversion) the 

angular correlation deviates from the Born approximation result if one takes full account 

of the nuclear charge number [^]. The angular correlation is then written as 

d^ri 1 d?7 , ^, 

' =---^ 1 + ecose . 2 
dii/dcosc) 2 aE 

Again the differential pair conversion coefficient ^ and the anisotropy coefficient e are 
computed numerically by employing the scattering solutions of the Dirac equation and 
assuming an extended nucleus. Using these Coulomb-distorted plane waves [^] internal 



pair conversion results in a strong dependence on the nuclear charge number. In Fig. 3 
the anisotropy coefficient is plotted versus the kinetic positron energy for the 1.76 MeV EO 
transition in Zr {Z = 40) (full curve) and U {Z = 92) (dotted curve). If the nuclear charge 
number is increased the angular correlation starts to deviate from the Born approximation 
result (dashed curve) and becomes more isotropic. 



3 Experimental setup 



As in the heavy-ion experiments |19, ^ 25|, we use two identical iron-free, orange-type 



(3- spectrometers facing each other with a common object point, at which a target 
wheel is placed (Fig. 4). Positrons (e^) emitted in the backward ('(9e+ = 110° — 145°) 
and electrons (e~) emitted in the forward hemisphere ("i^e- = 38° — 70°) are focussed 
onto the corresponding lepton detectors. Each lepton detector consists of an array of 
high-resolution Si PIN diodes, called PAGODA. Intrinsic to this setup is the capability 
of focussing, at a given field setting, only a certain momentum interval of e"^ or e~ by 
rejecting the opposite charge. This is of major advantage because the high f3~ background 
is suppressed completely on the e"*"-side, while the selection of only a certain band of 
e~-momenta on the e~-side enables operation with the strong sources required by the low 
branching ratios. 

The leptons are identified by matching their momentum, derived from the hit-point on 
the PAGODA and the spectrometer field setting, with their energy, measured by the PIN 
diodes. This efficiently suppresses backgrounds due to 7 rays. An additional coincidence 
requirement with the 511 keV annihilation radiation is thus not necessary for the e"*" 
identification. The momentum-energy matching in addition rejects the ~ 25% leptons 
backscattered from the PIN diodes almost completely, a unique feature of this setup. As 
mentioned above, for the lepton identification only events for which the energy-momentum 
relation is fulfilled are accepted, with the result that the remaining lepton misidentification 
is small and can be determined reliably. This is demonstrated in Fig. 5 for the case of 
the positron identification. As can be seen a clear signature for the focussed positrons is 
obtained. In such a spectrum all other particles, like electrons and 7 rays scattered from the 
vacuum vessel walls and the spectrometer coils as well as positrons backscattered from the 
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detectors, result in a broad continuous background distribution which can be determined 
quantitatively. The situation is very similar for the case of the electron identification. 

Each of the two identical lepton detection systems consists of 72 segmented, high- 
resolution, 1 mm thick. Si- PIN diodes which are mounted on 12 roofs, 6 detectors on 
each roof (see Fig. 4). Each detector chip has a trapezoidal shape (24 mm base, 16 
mm top, and 16 mm height) which is tilted by ~ 40° relative to the spectrometer axis. 
Each detector is also subdivided into three electrically separated sectors, each covering 
an azimuthal angular range of A0;ept = 20°. A matrix readout, each roof as well as 
each three neighboring sectors at the same azimuthal angle to one preamplifier, records 
information on energy and arrival time of a lepton hitting a detector. Cooling the detectors 
to ^ — 25°C, the e"'"e~-sum-energy and time resolution achieved amounts to ~15 keV and 
^4 ns (FWHM), respectively. 

The opening angle of the e"'"e~-pair, Og+e-; is measured directly within a range of 40° — 
180° in the laboratory. Using the 0-separation of the PAGODA's, this range can be sub- 
divided into ten angular bins (see table 1) with centroids in the range Oe+e- = 70°. . . 167°. 
These angular bins cover the same solid angle of 0.62 sr^, except the first and last, which 
cover only half of this value. This is due to simple combinatorics: for each of the angular 
bins all combinations between e"*" and e~ PAGODA columns with the same |A</)| are put 
together. For |A</)|=0° and 180° only one column on the other PAGODA contributes, but 
for e.g. |A</)|=20° both signs are possible and two columns of the other PAGODA con- 
tribute to this Ge+e--bin. The values given in table 1 have been calculated via a Monte 
Carlo simulation, which takes into account a realistic source spot size and the shadowing 
of the coils via ray tracing. The same code, but extended to handle the angular correlation 
of the emitted pairs, small-angle scattering and energy loss in the source, was also used 
for results shown later in this paper. 

Each pagoda array covers a maximum momentum acceptance of Ap/p =30% which 
corresponds to an energy interval of AE ~150 keV at a lepton energy of 300 keV. Within 
this momentum interval, the full-energy peak efficiency is 10% and 11% of Att, for electrons 
and positrons, respectively (see also Ref. p5[]). 



To cover an area in the E-^-E^ plane larger than given by the momentum acceptance 
of the spectrometers, the spectrometer currents are stepped up and down in a correlated 
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manner, where each step is done after the same hve-time interval. Thus an intrinsic 
correction for dead time effects due to different count rates is achieved. 

Under the assumption, that the ^-dependant longitudinal dispersion of the spectrome- 
ters cancels, the pair detection efficiency of the double-ORANGE setup can be written as 
follows: 

£ — Cspec{E-E, E/^, Igpec) X Canglei^e+e-) j (3) 

where Cgpec is the detection efficiency of the spectrometers for isotropic emission, depend- 
ing on energy distribution and spectrometer current sweep, and Cangie depends on the 
angular correlation of the leptons. However, for any given angular correlation, the relative 
number of counts in every angular bin of the setup does not depend on energy distribution 
or spectrometer fields. The angular correlation can thus be deduced from the relative 
intensities alone. The above mentioned assumption is justified, if leptons entering the 
spectrometers at any accepted angle 9 from the axis and passing the spectrometer will hit 
any of the pagoda's detectors. For a given, fixed spectrometer current this is true for the 
fiat top of the efficiency of ^ fa 30%. At the edges of the fiat top the efficiency drops fast 
to zero, and the efficiency in this case depends on 9. When the spectrometers are used 
in current sweep mode, the flat top region is extended over the whole acceptance range 
covered, and a 9 dependance is present only at the steep edges. 

4 Experimental results 

With the setup described above wc investigated the 1.76 MeV EO transition in ^°Zr and 
the 1.77 MeV Ml transition in ^°^Pb, using radioactive sources. 

4.1 The 1.76 MeV EO transition in ^^Zr 

The first excited 0+ state at 1.76 MeV of the 0+ nucleus ^°Zr is forbidden for normal 
7 decay and it's main decay channels are internal conversion, emitting electrons with 
sharp energy given by the transition energy minus the e~ binding energy 2/3) and 
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internal pair conversion 1/3)0) emitting e"'"e~-pairs with sharp sum energy given by 
the transition energy minus the equivalent of the pair's mass. The sum energy is thus 
Ey, = 1760keV— 1022 keV = 738 keV. This energy can be shared arbitrarily between e~ 
and e^, leading to a broad Ej\ distribution. Due to final state Coulomb interaction with 
the nucleus' field this distribution is shifted to higher e"*" energies and thus positive Ej\ 
values. Other decay modes of this transition have been studied in detail theoretically 
30|, |, |10|, HI as well as experimentally ^ 0, 



The 1.76 MeV transition decays to ~30% by IPC and is populated with a probability 
of 1.1 X 10~^ by the /3~ decay of ^°Ywith a Q- Value of 2.6 MeV. Since this isotope is 
rather short-lived (64 h), we used a ^°Sr source, where is produced in equilibrium by 
Strontium /5~ decay with a Q- value of 586 keV. Thus every e+e~ pair created via IPC decay 
of the first excited state of ^°Zr is accompanied by about 2 x 10^ electrons from /5~ decay. 
Especially the electrons from the ^°Sr decay cover exactly the range of investigated e"*" and 
e~ energies of :^150. . . 600 keV, determined by the low energy cutoff below 150 keV of our 
trigger electronics. To achieve satisfactory counting rates for a high statistics experiment, 
a very strong (in the order of 1 mCi) ^°Sr source had to be used, leading to total countrates 
on the e~ PAGODA of up to about 2 MHz as in the case of heavy-ion experiments. A very 
thin source had to be produced in order to keep broadening of the angular correlation by 
small-angle scattering of the outgoing leptons as small as possible. A source with 1.5yum 
Mylar foil plus a 0.2/im Aluminum (needed to keep the source at ground potential) was 
produced from commercially available ^°Sr solution^. 

Using this source we collected a number of 1.3x10^ e"'"e~-pairs in the sharp sum energy 
line of the IPC transition over a background of mainly chance coincidences of e+ from 
IPC and f3^ (Fig. 6a). However, the chance coincidence background can be determined 
quantitatively and subtraction yields the pure IPC spectrum (Fig. 6b). This spectrum 
is essentially background-free with a ratio of counts in the line (FWHM 15.5 keV) to the 
^ A very small part actually decays via multiple 7 emission. The strongest of these branches emits two 
7's back to back j2^. There is even a report of a very weak single-gamma branch made possible by 

coupling to the atomic shell electrons 

■^Our thanks to Dr. Briichlc and Co-workers, GST Nuclear Chemistry Division, for the production of 

this source 
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rest of the spectrum of 5.7:1. The tail on the low energy side is due to a small part of not 
suppressed backscattered leptons on either side. 

The high statistics in this experiment allows us to subdivide the data into 17 E/^ 
intervals of 66keV width and the ten opening angle bins of our spectrometer. Theory 
predicts for this EO transition an angular correlation of the form 

/(Ge+e-) =a- (l + eCOsGe+e-), (4) 

where a is given by the conversion coefficient, source strength, detection efficiency and 
duration of the measurement. To determine the angular correlation coefficient e the above 
given function is fitted to the measured data for each bin (Fig 7a), multiplying the 
first and last bin by two to correct for the smaller solid angle (rf. section However, 
the measured angular correlation includes also some small-angle scattering of the leptons 
in the source. To account for this the theoretical results were folded with small-angle 
scattering and the experimental response via a Monte Carlo simulation. The result of the 
simulation was fitted in the same manner as the measured data (Fig. 7b). Now the two 
values for e can be compared. The effective source thickness of 550 /ig/cm^ was calculated 
from the known thickness of the foils and with the assumption that the source material is 
pure. However, since the purity and the mean size of the crystallites of the active material 
are not well known, only a upper limit of about 1500 /ig/cm^ can be given. 

In Figure 8 we compare the anisotropy coefficients e determined from the measurement 
with the theoretical calculations adapted to our experimental conditions via Monte Carlo 
simulation for various effective source thicknesses indicated. Even for a 1500 yUg/cm^ thick 
source the theoretical values lie above the measured values except for very large \E/^\ values. 
This is a hint that the real anisotropy coefficient might be somewhat smaller (~9%) than 
the calculated values. However, within the uncertainties of the current analysis (dominated 
by the uncertainty of the source thickness) theory and experiment are consistent. The 
anisotropy given by older Born approximation results [with e{E/^ = 0) = 1] is inconsistent 
with the measurement. 
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4.2 The 1.77 MeV Ml transition in ^o^pb 

Using a rather weak (only a few yuCi) ^'^^Bi calibration source we were also able to measure 
e"'"e~-pairs from the weak IPC branch of the 1.77 MeV Ml transition in ^°^Pb, leading to 
a narrow e"'"e~-sum-energy line at E^, = 755 keV. Since we detected here only a few pairs 
per hour, the statistics of this measurement does not allow to subdivide into several E/^ 
bins. Due to the low rates practically no chance coincidences occur in this measurement, 
and the measured e"'"e~ sum energy spectrum shown in Fig. 9 is again of the same quality 
as that of Fig. 6b. 

However, even the low statistics of this measurement allows us to subdivide into the 
ten opening angle bins and determine the angular correlation of the pairs. The theoretical 
prediction for /(Gg+e-) is in this case more complex, but the form given in eq. ^ is the first 
order in a cos(0e+e-) expansion of the correct result, so again the same fitting procedure 
can be used. As can be seen in Fig. 10, the situation is very different from the EO case 
reported in the last section: while for the EO case the intensity decreases for larger Og+e- 
values, in this case the intensity increases in the Qe+e- range covered by our setup, thus 
yielding a negative value of e = —0.12 ± 0.06 for the anisotropy coefficient. This is in 
strong disagreement with the theoretical results derived using the Born approximation, 
but in good agreement with the most recent theoretical results (see Fig. 1). 

These new results have some impact on estimates for the expected strength of the 
accompanying 7 line, if some of the previously reported sum energy lines found in heavy 
ion experiments could be explained by IPC processes. According to these findings the pair 
efficiency of our setup is much higher for angular correlations having the pattern of the Ml 
case than for angular correlations derived from Born approximation. This is particularly 
true for lines appearing only in the 180° bin of the setup, where the Doppler shifts of the 
two leptons almost cancel: For this specific situation the detection efficiency is an order 
of magnitude larger for e = —0.12 than for e = 1. We would thus expect a factor of ten 
weaker 7 line for e = —0.12 than in the old estimate p4|. Only a detailed re-analysis of 



the old data p4[ can clarify whether IPC could be the origin of these weak lines. It should 



*No ^"^Bi source material was commercially available to produce a much stronger source in the same 
way as with ^"^Sr 
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be noted in this context that recent results, obtained from high-resolution Doppler-shift 
spectroscopy in heavy-ion collisions at the Coulomb barrier, have indeed revealed that 
weak e+e~-sum-energy lines from IPC transitions of moving emitters can appear in the 



measured spectra p4|. 



5 Summary 

The measurements presented in this work prove the capability of our setup to detect the 
sharp e"'"e~ sum energy line of the IPC decay of excited heavy nuclei despite the several 
orders of magnitude more intense background of electrons and 7 rays from the same source. 
The measurement with the ^°Sr source clearly shows that this is even the case for very 
high count rates in the MHz range on the electron detectors, as they are typical also for 
in-beam measurements (see e.g. Ref. p5| ). 

The results of the first energy and angle resolved measurements of IPC pairs — made 
more than 60 years after the first theoretical prediction of this decay branch — prove that 
Born approximation is not suitable for the treatment of IPC transitions of only a few times 
the threshold in heavy nuclei even for the most simple EO case and by no means for mag- 
netic multipole transitions. The theoretical results for the angular correlation presented 
here, taking into account all relativistic corrections, the finite size of the nucleus and also 
the magnetic substates for magnetic transitions, are in agreement with the experimental 
results shown in this paper. The combined theoretical and experimental efforts of our 
groups have thus led to a better understanding of a fundamental decay channel of excited 
nuclei that was assumed to be very well understood for decades. 
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TABLES 



Table 1: e+e opening angle bins of the ORANGE setup 
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FIGURE CAPTIONS 



Fig.l: Doubly differential pair conversion coefficient plotted versus the opening angle 
G of the electron-positron pair for an Ml transition in ^°^Pb. The transition energy 
amounts to 1.77 MeV, the kinetic positron energy was fixed to 450 keV. The dotted line 
represents the Born approximation result, which is valid for small nuclear charge numbers. 
In contrast to this, the DWBA calculation yields a nearly isotropic distribution. One can 
also deduce the influence of the nuclear size (full line) which is compared to the point 
nucleus result (dashed line). 

Fig. 2: Angular correlation for a Ml transition with transition energy 1.77 MeV. The 
positron energy is fixed to 370 keV (symmetric splitting of the nuclear transition energy 
on electron and positron). In order to study the effect of the nuclear charge the angular 
distribution is plotted for the nuclear charge numbers Z — 92 (full curve) , Z — 82 (dashed 
curve), Z — AO (dash-dotted line), as well as for the Born approximation result (dotted 
line) . 

Fig. 3: Anisotropy coefficient e plotted versus the energy difference E/^ for EO conver- 
sion assuming a transition energy of 1.76 MeV and three nuclear charge numbers Z = 
(Born approximation result, dashed line), Z = AO (full line) and Z = 92 (dotted line). The 
Z — 40 curve is the prediction of the anisotropy coefficient for the EO transition in^^Zr. 

Fig. 4: Schematic view of the new double-Orange setup. Each of the /9-spectrometers 
is equipped with a /9-multidetector system of 72 Si (PIN) diodes (e^-Pagodas) . The 
forward spectrometer is surrounded by 18 position-sensitive heavy- ion detectors (PPAC), 
and contains a further PPAC detector in its center. Also shown is the rotating target 

wheel and the Gc(i) 7-ray detector. The only change made for the current experiment is 
that the target wheel was used in a stationary mode, on a target position of which the 
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radioactive sources were mounted. 

Fig.5. Experimental signature of the detected positrons in the ^°Zr measurement. The 
abscissa shows the normahzed difference between the momentum of the focused positrons 
calculated from the pulse height of the Si detectors and their momentum determined from 
the magnetic field of the spectrometer. Only events falling in a narrow window centered 
around zero are accepted in the analysis. The background is dominated by electrons from 
the five orders of magnitude more intense (3~ rays, reaching the positron detectors after 
multiple scattering in the spectrometer. 

Fig.6: e"'"e~-sum-energy spectrum from the ^°Sr source. Left part: raw spectrum, 
overlayed with the chance coincidence background due to the high (1. . . 2MHz) /3~ count 
rate determined from not prompt coincidences. Right part: Sum energy spectrum after 
subtraction of chance coincidences. The small low energy tail is due to a small fraction of 
backscattered leptons not rejected by the lepton identification (see text). The sum energy 
line from the 1.76MeV IPC decay in ^°Zr is measured essentially background- free despite 
the five orders of magnitude higher /3~ background from the decay of ^'^Sr and ^°Zr. The 
resolution is 15.5 keV (FWHM), the total content of the spectrum 1.59x10^ 

Fig. 7: An illustrative example of the determination of the angular correlation coef- 
ficient e for both, measurement (left part) and theory (right part), which are treated in 
the same manner: The function / = a(l -|- ecosGe+e ) is fitted to the number of counts 
in the angular bins of the ORANGE setup at an energy difference centered at zero. The 
theoretical result is gained by folding the theoretical prediction with small-angle scattering 
in the source and the acceptance of the ORANGE setup by a Monte Garlo simulation. 
The error quoted in the experimental value of e = 0.83 ±0.03 is solely of statistical origin. 

Fig. 8: The extracted angular correlation coefficient e for IPG decay of the 1.76 MeV 
transition in ^°Zr as a function of the energy difference £'a of the emitted e+e~-pairs. The 
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theoretical prediction (solid line) was folded with small-angle scattering in a source with 
a thickness of 750/ig/cm^ (dashed hne) and 1500 //g/cm^ (dotted hne). The points show 
the experimental result. 

Fig.9: e+e~ sum energy spectrum measured with a weak ^^'''Bi source. Due to the low 
count rates, chance coincidences of leptons from IC and P'^ decay branches are negligible. 
Again the sharp sum energy peak from the IPC decay of the 1.77 MeV Ml transition in 
^°^Pb is measured essentially background-free, even though a huge (five orders of magni- 
tude larger) background of 7's, e~ from IC and /3+ is present. 

Fig. 10: Measured angular distribution of the IPC line from the 1.77 MeV Ml transition 
in ^°^Pb. Again the intensity in the first and last bin have been multiplied by two to account 
for their smaller solid angle. A fit done in the same manner as in Fig. 7 yields a negative 
anisotropy coefficient e = —0.12 ± 0.06, i.e. the minimum of the intensity is at 90° rather 
than 180°. 
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